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Microbial colonies are experimental model systems for studying the colonization of new territory 
by biological species through range expansion. We study a generalization of the two-species Eden 
model, which incorporates local frequency-dependent selection, in order to analyze how social inter¬ 
actions between two species influence surface roughness of growing microbial colonies. The model 
includes several classical scenarios from game theory. We then concentrate on an expanding public 
goods game, where either cooperators or defectors take over the front depending on the system pa¬ 
rameters. We analyze in detail the critical behavior of the nonequilibrium phase transition between 
global cooperation and defection and thereby identify a new universality class of phase transitions 
dealing with absorbing states. At the transition, the number of boundaries separating sectors decays 
with a novel power law in time and their superdiffusive motion crosses over from Eden scaling to a 
nearly ballistic regime. In parallel, the width of the front initially obeys Eden roughening and, at 
later times, passes over to selective roughening. 

PACS numbers: 02.50.Le, 68.35.Ct, 68.35.Rh, 87.18.Hf, 87.23.-n 


I. INTRODUCTION 

Living species are usually confined to their territory, a 
spatial region defined by geographical borders, climate, 
or other environmental constraints. Uninhabited regions 
are colonized through range expansion, where individuals 
reproduce and disperse at the front of their territory [T]. 
This process is seen in biological invasions [2] , as a result 
of shifting climate zones mg, during colonizations in 
our own species’ history EHH], tumor growth [nisiiin], 
and biofilm growth mHig. Evidently, expansions occur 
on very different spatial (micrometers to 10^ meters) and 
temporal (hours to millennia) scales. 

In this article we aim to characterize range expansion 
under the influence of short-range “social interactions” 
of individuals at the front. Such interactions are present 
if success in reproduction depends on the presence of 
nearby individuals of the own and/or other species. Here, 
we set up a model for the expanding front based on evo¬ 
lutionary game theory [141116) and investigate its rough¬ 
ening dynamics for two interacting species. Besides ex¬ 
ploring an interesting non-equilibrium growth process, we 
hope to contribute to interpreting experiments on range 
expansion in multi-species colonies of simple organisms. 

In experiments, microbial growth is excellently suited 
to study range expansion and other processes in pop¬ 
ulation dynamics and evolution such as spatial spread 
of infections and adaptation to an environment (see for 
example Ref. [HI)- Microbes reproduce fast, their en¬ 
vironment and genotype can be controlled, and exper¬ 
imental conditions are easily reproducible. Grown in a 
Petri dish, the spatial patterns of single-species microbial 
colonies have long been a rich field of study . The 

observed patterns crucially depend on motility, availabil¬ 
ity of nutrients, and the growth medium, to name but a 
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few. However, even under conditions of negligible motil- 
ity and abundant nutrients a colony’s front is rough and 
has interesting statistical properties [23H25] . 

Multi-species colonies are composed of more than one 
species and show additional intriguing features, even if 
the species are identical except for a marker |26) . During 
reproduction they keep their marker but compete with 
other species for space at the front. Thereby, sectors of 
single species form, which are separated by boundaries. 
Their statistical and dynamic properties are determined 
by the evolving roughness of the expanding front EZl un¬ 
usually, when two or more species of microbes live in 
a common environment, they influence each other dur¬ 
ing reproduction. In particular, reproductive success of 
any species, also called its fitness, depends on the pop¬ 
ulation sizes of all the species. This constitutes “social 
interactions” between the species commonly referred to 
as frequency-dependent selection. Research in the field 
has initiated a wealth of fascinating experiments uni Hi- 
Ill] either in well-mixed liquid culture without any spatial 
order [31] [36] or in the Petri dish, where the populations 
are spatially structured [IZlIlHj. Many of the experimen¬ 
tal observations can be discussed within the framework 
of evolutionary game theory dMIS]. For example, light 
has been shed on a long-standing theoretical question in 
evolution [H HU [Ml SO]: Why do individuals cooperate 
if non-cooperators can exploit them? Literature empha¬ 
sizes the importance of a population to be structured 
in groups [HI 1411145] . for example, by spatial distance. 
While within a single group cooperators are always in¬ 
ferior to non-cooperating defectors, if the latter interact 
only with their neighborhood, distant large groups of co- 
operators will ultimately outcompete defectors. Some 
models also stress the central role of demographic fluctu¬ 
ations and of populations growing in size [US?]. Both, 
experiments and theory, explain the advantage of coop¬ 
erators during colony growth by their ability to locally 
advance faster [Ml EHl SS], vividly termed “survival of 
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the fastest” [55] , 

Cooperation between between nearby cells is often me¬ 
diated by some biochemical compound (a public good) 
which the microbes release into the extracellular envi¬ 
ronment. This compound then promotes reproduction of 
neighboring cells. In general, a released substance may 
act beneficial or detrimental to other individuals, also 
depending on their species, and implies some cost to the 
producer. Examples include secretion of digestive inver- 
tase to break down sucrose [Ml [5DHM] , siderophores to 
scavenge iron from the environment [511 m m] , poly¬ 
mers which support biofilms UU 1551IM] . release of tox¬ 
ins [SllllH] (sometimes through lysis [591160]), surfactants 
which facilitate swarming m, and the exchange of amino 
acids [57] . 

This plethora of biochemical compounds, released by 
cells and affecting nearby cells, implies a wealth of spe¬ 
cific features, which certainly are not covered by a sin¬ 
gle model. However, since the released biomolecules 
usually mediate short-range interactions between indi¬ 
viduals, properties on large scales should be indepen¬ 
dent of microscopic details. Hence, we formulate a sim¬ 
ple model which captures the essence of an interaction 
while ignoring complicated details. The classical Eden 
model [HU [55], a simple growth process on a lattice, has 
been used successfully to mimic growing cell colonies.. It 
generates a cluster (the colony), the surface of which ex¬ 
hibits scaling properties also found for expanding fronts 
of microbial colonies [53] [55] . Extended to two identi¬ 
cally growing but still distinguishable species, it gener¬ 
ates sectors occupied by a single species only [28l 155] . 
Indeed, this behavior is found for two-species micro¬ 
bial colonies [55]. Moreover, boundaries between sectors 
move superdiffusively as in the experiments. 

In this article we explore a generalization of the two- 
species Eden model, which incorporates local frequency- 
dependent selection. We thereby aim to analyze how 
social interactions influence surface roughness of grow¬ 
ing microbial colonies. We set up an expanding public 
goods game, where either cooperators or defectors take 
over the front depending on the system parameters [M]- 
[15]. Right at the transition the front displays critical 
behavior, which we analyze in detail. In particular, we 
establish that our model belongs to a new universality 
class of phase transitions dealing with absorbing states. 
At the transition, the number of boundaries separating 
sectors decays with a novel power law in time and their 
superdiffusive motion crosses over from Eden scaling to 
a nearly ballistic regime. In parallel, the width of the 
front initially obeys Eden roughening and, at later times, 
passes over to what we call selective roughening. 

The remainder of this article is organized as follows. 
To analyze multi-species microbial colony growth, we in¬ 
troduce the Eden model with frequency-dependent selec¬ 
tion in Sec. jH] and analyze its phenomenology in Sec. |HI| 
We then concentrate on the expanding public goods game 
with its social dilemma in Sec. |IV| and analyze the critical 
behavior at the transition between long-term cooperation 
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FIG. 1. Two-species Eden model with frequency-dependent 
selection on a hexagonal lattice. The bacterial colony grows 
from the bottom line (lattice sites with narrow black edge) 
of length L, where individuals of species 1 (blue) and species 
2 (red) occupy the lattice sites. The colony expands into 
the empty, infinitely-extended half-space. Individuals capable 
of reproduction (indicated by bold edges), have at least one 
empty lattice site as a nearest neighbor. In a reproduction 
event, one of these empty neighboring sites {i,j) is chosen 
with equal probability and the reproducing individual changes 
the corresponding state Sij to its own state 1 or 2. Each 
individual has its own reproduction rate bij given in Eq. 0 . 
For example, the reproduction rate of the individual at site 
(3,3) (bold black edges) is 63,3 = 62 + 3T -I- 21^. Along the 
transverse direction periodic boundary conditions apply. 


and long-term defection by applying statistical analysis. 
Finally, we discuss and summarize our findings in Sec.jV] 


II. EDEN MODEL WITH 
EREQUENCY-DEPENDENT SELECTION 

In this work we employ a lattice model (see Fig. 
to analyze range expansion at rough fronts under the in¬ 
fluence of frequency-dependent selection. We set up a 
cellular automaton on a two-dimensional hexagonaQ lat¬ 
tice of transverse extension L and an infinite longitudinal 
extension. Periodic boundary conditions are applied in 
the transverse direction. The state {s} of the system at 
time t is specified by the state variables Sij of lattice sites 
(f,j). Consider a system with two species (extension to 
more species is straightforward). For any time t, any site 
(j,j) is either empty {sij = 0) or occupied by an indi¬ 
vidual of either species 1 {sij = 1) or 2 {sij = 2). All 
individuals which have at least one free nearest neighbor 
site can reproduce. To perform a reproduction step, we 
choose one of these fecund individuals with a probability 
proportional to its reproduction rate (see below) and a 


^ On a square lattice it is impossible to enclose a cluster A within 
a cluster B, which only contains nearest neighbor sites of cluster 
A. On a hexagonal lattice this is possible. 
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new individual of the same species is placed with equal 
probability on one of the free neighboring sites. 

In contrast to the Eden model [S2] and some of its two- 
species generalizations [28l [67] , reproduction rates in our 
model depend on the states of the nearest-neighbor sites. 
Let ni and n 2 denote the number of nearest neighbors of 
species 1 and 2 , respectively, then the reproduction rate 
of an individual at lattice site (*,j) is 

{ 0 if Sij has no free neighbors, 

6 ° -I- niTZ -I- n2S if Sij = 1, (1) 

62 + niT + n2'P if Si j = 2 . 

Here, b\ and 63 are the respective contributions to the 

reproduction rates of species 1 and 2 , which are indepen¬ 
dent of the states of their nearest neighbors. Frequency- 
dependent selection is introduced through the parameters 
TZ, S, T, and V. 

With the reproduction rates bij we implement a ran¬ 
dom sequential update of the system using a simplified 
version of the Gillespie algorithm [ 68 ]. The overall re¬ 
production rate of the population is btot ■= j ^i,j 
an individual at site (z,j) is selected to reproduce with 
probability bij/btot- We then choose one of the empty 
nearest-neighbor sites of the reproducing individual at 
random and place there a new individual of the same 
species. This implies that there are no mutations. Since 
the mean time until the next reproduction event is 6 ^^, 
we update time hy t ^ t + b^J^ after each reproduction 
event. We assume that individuals do not die and that 
they are immobile. Therefore, any site with Sij ^ 0 
remains in its specific state indefinitely. As initial condi¬ 
tion we occupy all sites of an initial line randomly, but in 
equal parts, with species 1 and 2 , if not stated otherwise. 

The formulated model generalizes version C of the 
Eden model, introduced by Jullien and Botet [M], to a 
two-species system. We already applied a similar model 
to range expansion without frequency-dependent selec¬ 
tion but included the possibility of mutations ]nZ]- If 
S, T, and V are zero, our model reduces to that of Saito 
and Miiller-Krumbhaar [28] , however they used a square 
lattice. Since diffusion is not included in the model, con¬ 
figurations and patterns behind the front are frozen. This 
corresponds to observations in microbial experiments on 
range expansion [26l [60] . 

In game theory the parameters TZ, S, T, and V 
from Eq. Q define the payoff matrix of a two-strategy 
game [Tllll6| . Different scenarios, some well known in 
game theory, are implemented if we set these parameters 
accordingly. 

III. PHENOMENOLOGY 

We now describe some generic examples of our model 
for growing microbial colonies (see Fig. [^, which emerge 
for typical parameter settings, and discuss their charac¬ 
teristic features. We then concentrate on so-called so- 



FIG. 2. Growth patterns of our model for different param¬ 
eters, which correspond to typical settings. Time is always 
t = 15 and lattice size is L = 200. Gooperators are depicted in 
blue and defectors in red. If not stated otherwise 62 = ^2 = I, 
TZ = V — S = T =0, and the initial ratio of both species is 
1 : 1 . (a) neutral growth, (b) selective advantage for species 2 
(62 = 1.5), which occupies 10% of initial sites, (c) coordina¬ 
tion game {TZ = P = 1), (d) snowdrift game (5 = T = 1), 
and (e) public goods game {TZ = 0.1,T = 1.1). 


cial dilemmas, where one species (defectors) exploits the 
other one (cooperators). Due to the spatial extent of our 
system, cooperators are able to outcompete defectors in 
a defined parameter region. This is in contrast to a single 
group, where all members interact with each other and, 
therefore, cooperators are always inferior to defectors. 

In the simplest case selection is frequency-independent, 
TZ — S = T = 1^ — 0, and both species reproduce with 
the same rate, = &q [see Fig. [^a)]. For this “neutral” 
setting we observe roughening of the front typical for the 
Eden model [H2E3]. Simultaneously, sectors composed of 
a single species merge and thereby coarsen [28l [ 66 ] ■ This 
inherent process happens according to the following sce¬ 
nario. If the tips of two advancing boundary lines meet, 
they annihilateand the enclosed sector loses contact to 
the front. Consequently, the number of boundaries and 
sectors can only decrease. Sectors repeatedly coarsen as 
they merge in these events. When all boundaries have 
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vanished, the front “has fixed” to a single species keeps 
on expanding. In finite systems, L < oo, fixation to a 
single species always occurs since in our stochastic model 
there is always a finite rate at which boundaries annihi¬ 
late. Hence, two absorbing states exist. Eventually, the 
expanding front will fix either to species 1 or species 2 . 

In Fig. [^b), species 2 has a larger reproduction rate, 
6 q > and therefore a constant selective advantage. 
As the front expands faster at locations where it is com¬ 
posed of species 2 , the roughness of the front increases. 
Indentations of the front usually are caused by sectors of 
species 1, whereas species 2 creates bulges. Furthermore, 
boundaries are biased such that sectors of species 2 widen 
while sectors of species 1 shrink laterally. Hence, sectors 
of species 2 merge and coarsen quickly. Eventually, the 
expanding front will fix to species 2 , which has almost 
happened in Fig. [^b). 

If the reproduction rates depend on the state of near¬ 
est neighbors (frequency-dependent selection), new pat¬ 
terns arise. In this article we are mainly interested in 
frequency-dependent selection and therefore set = 
62 = 1 from here on. We now discuss some interesting 
cases, see Fig. |^c)-(d), which correspond to well-known 
settings in game theory [I4HI6] . 

In coordination games (72. > T and V > S), see 
Fig-i: c), the front expands slower near boundaries than 
in the centers of sectors. Therefore, indentations in the 
front are typically found where boundaries currently are 
or recently have been. After sectors have coarsened for 
some time, most individuals are located inside sectors. 
Therefore, most of them only have neighbors of their own 
kind, which raises the average reproduction rate and the 
overall front advances faster. 

In snowdrift games (72 < T and V < S), see Fig.j^d), 
the front expands faster near boundaries. They annihi¬ 
late less frequent as compared to Fig. [^c) since narrow 
sectors grow faster. Boundaries are also strongly twisted 
and associated with bulges of the front. 

In this article we are mostly interested in social dilem¬ 
mas where one species (called cooperators) raises the 
reproduction rate of all neighbors regardless of their 
species. The increased reproduction rate is called a pub¬ 
lic good in game theory, since it is of benefit to all nearby 
individuals, but it also costs resources. In contrast, the 
other species (called defectors) takes advantage of the 
public good for its own reproduction, but does not con¬ 
tribute to the reproduction of its neighbors in the same 
way. Defectors save resources for their own reproduction 
and therefore have an advantage. In this scenario de¬ 
fectors do not at all contribute to reproduction of their 
neighbors and, therefore, we set 5 = 7^ = 0 and just 
vary 72 and T- According to Eq. Q, these two parame¬ 
ters increase the respective reproduction rates of coopera¬ 
tors and defectors if they have cooperating neighbors. In 
Fig-i: e) we present a setting, where species 2 (defectors) 
rapidly takes over large parts of the front. Defectors ben¬ 
efit from the initially large number of boundaries, where 
they take advantage of nearby cooperators, and conquer 



FIG. 3. Schematics of possible scenarios in an expanding 
public goods game. Depending on the parameters 72, 5, T, 
and V of Eq. 0 . cooperators (C, blue) and defectors {D, red) 
advance with different speeds, as indicated by arrows, (a) 
Cooperators outrun the trailing defectors. From an advanced 
position along the front cooperators can then expand laterally 
and take over the front, (b) A thin layer of defectors keeps 
up with the cooperators’ sector and, eventually, completely 
covers the cooperators. 


most of the front. Only cooperators, living in sufficiently 
large sectors, can keep up with the front during this early 
period and may then take over the front, depending on 
the parameter values. 

In a situation like this, it is not a priori clear if the 
front eventually fixes either to cooperators or to defec¬ 
tors. Depending on the values of 72 and T, cooperators 
can either outrun defectors and, from their advanced po¬ 
sition at the front, overgrow their competitors laterally, 
see Fig.j^a). Or, defectors cover cooperators with a thin 
layer and thereby take over the front, see Fig.j^b). Close 
to the transition between both scenarios, the front dis¬ 
plays increasing roughness since both species are able to 
take over while their fronts grow with different speeds. To 
characterize this transition quantitatively, we performed 
extensive simulations and applied methods from surface 
roughening [631 EDHig and the theory of phase transi¬ 
tions dealing with absorbing states [71 Eg. 


IV. EXPANDING PUBLIC GOODS GAME: 

CRITICAL BEHAVIOR 

In this section we quantitatively analyze the transi¬ 
tion between long-term cooperation and long-term defec¬ 
tion for an expanding public goods game. As the transi¬ 
tion is approached, several observables show critical scal¬ 
ing EHESl- Following our earlier work EZ], we perform 
finite-size scaling to localize the transition. Furthermore, 
we determine critical exponents and thereby establish a 
new universality class for the transition between the two 
adsorbing states. In the vicinity of the transition we also 
study the dynamics of the sector boundaries including 
the decline of their mean number during coarsening and 
their superdiffusive motion as well as the roughening of 
the expanding front. 
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FIG. 4. Probability Pax of the front to fix to cooperators 
plotted versus T for several system sizes L at TZ = 0.1. Er¬ 
ror bars give the standard error of the mean for each data 
point. Lines are fits of the data points to the Fermi function, 
1 /exp(c(T — T 1 / 2 ) + 1), where 71/2 and c are fit parameters. 
Inset: The transition point 7i/2(L) relative to the fitted crit¬ 
ical value Tc ~ 1.58 (black crosses) follows a power law in 
L —>■ 00 : Tc — 71/2(7/) = (j 4/L)^^‘*'^ (black line). 


A. Finite-size scaling and phase diagram 


FIG. 5. Phase diagram of the two-species public goods game 
with range expansion. Parameter regimes, where the expand¬ 
ing front fixes either to cooperators or defectors in large sys¬ 
tems with L —>■ 00 , are indicated by blue and red shade, re¬ 
spectively. Red crosses are from simulations with system size 
L = 1000, where defectors always outcompeted cooperators 
(Pfix = 0) while blue dots identify events where cooperators 
survived (Pax > 0). The black diamonds indicate critical 
points Tc for P —> 00 determined from finite-size scaling (see 
Eq. § and Fig. 1^. Note that in finite systems defectors have 
an advantage in a larger parameter region. 


As boundaries merge, sectors coarsen and the system 
progresses towards one of the two absorbing states. At 
finite system sizes L this is a stochastic process and both 
adsorbing states are reached with a certain probability. 
However, in the thermodynamic limit, L —>■ 00 , the mag¬ 
nitude of fluctuations relative to the mean value goes to 
zero and one of the absorbing states is reached with cer¬ 
tainty. We now use the method of finite-size scaling to 
determine the transition point between both states Ca¬ 
in Figure we present the probability Pfix that the 
front fixes to cooperators and plot it versus T for several 
lattice sizes L at 7^ = 0.1. We distinguish two regimes: 
one where cooperators dominate (Pgx(T, A) > 5 ) and 
one where defectors take over. We locate the transition 
point atTi/ 2 (A) by Pflx(Ti/ 2 (A), A) = ^. As L ^ 00 , Pfix 
converges to a step function, since in infinite systems the 
absorbing states are reached with certainty. The step is 
positioned at Tc '■= limi_).oo 7 i/ 2 (A). From the theory 
of critical scaling applied to absorbing states, we expect 
that close to the critical point Tc the states of the lattice 
sites are correlated on the transverse distance f_L. Ap¬ 
proaching Tc, C-L diverges as where A ■= Tc — T 

is the distance to the critical point and is a critical ex¬ 
ponent. For finite systems, an absorbing state is reached 
if 

A « a ~ lAf"" . (2) 

The transition occurs at 71/2 (A) and rearranging Eq. ([^, 


we obtain 

ri/2(A) « Tc - (A/A)!/-^- . (3) 

The characteristic length A is related to the microscopic 
length scale, which here is the lattice constant, and de¬ 
tails of our model. It is not important to the following 
analysis. The inset of Fig.|^shows the best fit of our data 
to Eq. which yields the critical exponent « 4.2 
and the critical point Tc ~ 1.58 at 77 = 0.1. 

The above procedure can be repeated for different val¬ 
ues of 77 to map out the phase diagram (see Fig. [^. 
One realizes that the benefit of cooperators from their 
own species, 77, has a much more pronounced influ¬ 
ence on the final state than the defectors’ benefit from 
cooperating neighbors, T- This makes sense, since at 
large times 7^1 the front contains large single-species 
sectors. Hence, the number of sector boundaries Nb, 
where defectors can benefit from cooperators, is small: 
Nb Ti A. Therefore, almost all cooperators have coop¬ 
erating neighbors, while only a few defectors have this 
advantage. This is an example of “preferential assort¬ 
ment” , where the benefit of cooperation is almost entirely 
available to other cooperators HU [311131 msii El US] . 
So, for a wide range of parameter combinations coopera¬ 
tors can indeed outcompete defectors. However, for large 
enough T the dynamics at the boundaries still determines 
the final state of the front and defectors outcompete co- 
operators. 














6 



FIG. 6. Mean time to fixation as a function of defector 
benefit T for various system sizes L at 7?. = 0.1. The inset 
depicts the rescaled fixation time, where A = 7^ — 7i All 
data collapse on a single master curve for critical exponents 
iz|| = 3.5 ± 0.1 and v±_ = 4.2 ± 0.1 and critical point Tc = 
1.59 ± 0.03. Accordingly, t^x grows with L and the position 
of its maximum approaches Tc for L —> oo. 


B. Critical exponents of the phase transition 



FIG. 7. Survival probability Pc of cooperators starting from 
a single site plotted versus time for several values of T. Other 
parameters are L = 1000 and TZ = 0.1. For T > Tc, Pc 
decreases exponentially, while for T < T the front fixes to 
cooperators with a non-zero probability. At the transition 
point Tc, the survival probability decays in time with a power 
law with exponent P'/v\\ = 2.2 ± 0.1. Inset: The number of 
cooperator sites at the front, Nc, decays exponentially in time 
for T > Tc and is non-monotonic for T < T- At the transition 
Nc decreases with a power law with critical exponent 0 = 
1.3 ±0.1. 


In the previous section IV A we already encountered 
the critical exponent iy±. We now continue to determine 
further critical exponents of the phase transition. These 
exponents are universal. They only depend on the dimen¬ 
sion of the system, the number of components of the order 
parameter, and symmetries of the model [71173]. They 
are independent of microscopic details and do not vary 
along a phase transition line. Our system has properties 
similar to “compact directed percolation” (CDP) [76] . 
This is a stochastic process with a flat front, which also 
has two distinct absorbing states. Using this similarity, 
we proceed by determining critical exponents, which are 
known for CDP m, and compare both models. 


At the critical transition, T = Tc, none of the two 
species has an advantage. Heterogeneous fronts, com¬ 
posed of more than one sector, exist for long times be¬ 
fore the front fixes to one of the absorbing single-species 
states. This can be quantified by the mean time to fix¬ 
ation, tfix, presented in Fig. [^ The data show that the 
fixation time has a maximum, the position and value of 
which grow with system size. 


Along the longitudinal direction, in which the front 
propagates, states are correlated on the longitudinal dis¬ 
tance ^ 11 . As before, we expect it to scale like | close 
to the transition. Since the front propagates with a mean 
velocity, is proportional to a correlation time. Close to 
Tc this time becomes very long, which is known as critical 
slowing down. Substituting Eq. § into ^11 ~ |A| , we 


find the scaling relation 

(4) 

We expect the mean fixation time to be proportional to 
the correlation time ^ Therefore, in the inset of 
Fig. y we plot ffix rescaled by versus A rescaled by 

_ All curves of the main plot collapse on a single 
master curve for Tc = 1.59 ±0.03, and critical exponents 
i/_L = 4.2 ± 0.1 and z/y = 3.5 ± 0.1. The values of Tc and 
i/_L are in good agreement with our fit to Eq. So, 
fixation of the front is determined by the characteristic 
time 


Tfix ~ Cll ~ . (5) 

Two more critical exponents right at the transition are 
related to the survival probability of one species or state, 
which initially occupies a single site while all the other 
sites are occupied by the other state. We choose a single 
cooperator site in a line of defectors and determine the 
probability Pc (t) that after time t there are still cooper¬ 
ators at the front and also calculate the average number 
of cooperator sites at the front, Nc{t). In Fig. [^we plot 
both quantities versus time for different defector bene¬ 
fit T. At the transition situated between T = 1.55 and 
1 .6, we find that both Pq and Nq (see inset) decay with 
power laws in time: Pcit) and Nc{t) ~ . 

The respective best fits yield P'/i^\\ = 2.2 ± 0.1 and 0 = 
1.3 ± 0.1. Using our result for zzy, we find (3' = 7.7 ± 0.6. 
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TABLE I. Critical exponents for the phase transition to the 
absorbing states (either long-term global defection or long¬ 
term global cooperation) for the expanding public goods game 
with frequency dependent selection. 



l/|| 

P 

P' 

e 

4.2 ±0.1 

3.5 ±0.1 

0 

7.7 ±0.6 

1.3 ±0.1 


In general, in phase transitions to absorbing states the 
critical exponent /3 governs the stationary density of “ac¬ 
tive sites”, when approaching the transition | 73 j . In our 
case, the “active” sites can either be cooperators or defec¬ 
tors. Since the stationary state is either an all cooperator 
or an all defector front, the density ~ |7)) — T\^ jumps 
from 0 to 1 and, hence, j3 is 0. Our results for all the 
critical exponents are summarized in Table |l] 

The set of critical exponents determines the univer¬ 
sality class of a phase transition. To our knowledge no 
other non-equilibrium transition has been found, which 
shares the same set of exponents. Hence, the transition 
between long-term cooperation and long-term defection 
in our expanding public goods game with frequency de¬ 
pendent selection constitutes a new universality class. 

C. Dynamics of boundaries 

In this section we investigate the dynamics of the 
boundaries which separates sectors of cooperators and 
defectors from each other. In rough fronts the local front 
orientation is tilted against the main growth direction. 
When the front grows further, this tilt directs the move¬ 
ment of boundaries [5^ [55]. Ultimately, when two of 
them meet, they annihilate. In Fig. we plot their mean 
number Ni, versus time for several values of the defector 
benefit T. Right at the transition (dashed black line), 
shows a power law decay. We now discuss the different 
regimes in Fig. 

For neutral systems, where frequency-dependent selec¬ 
tion is absent (T = 77. = 0 ), any inclination of the front is 
created by stochastic surface or Eden roughening [ 5511101 “ 
[ 72 ]. The surface undulations obey KPZ-scaling m and 
thereby drive the decay of iVf, [551HH] ■ Boundaries move 
superdiffusively along the front with a mean-square dis¬ 
placement proportional to [57]. On average, they 

annihilate after having traveled the mean distance L/Nj, 
between the boundaries, for which they need the time 
^ (L/7Vb)^/^. Hence, boundaries annihilate with a rate 
proportional to , which implies 

Nb ~ . ( 6 ) 

So, the number of boundaries decreases as 

Nb{t) ~ r 2/3 , ( 7 ) 

as already observed by Saito and Miiller-Krumbhaar [55] . 
This power law is excellently reproduced by our simu- 



FIG. 8. Mean number of sector boundaries Ni, plotted as 
a function of time for several values of 77 Other parameters 
are L = 1000 and 77 = 0.1. Solid black line: Ni,{t) ~ 
as predicted in Ref. [28] for neutral growth (T = 77 = 0). 
At large times the simulated curve deviates from the power 
law due to finite system size. Dashed black line: Close to 
the critical point T = 77 ~ 1 .6, we find a power-law decay 

lations in the case of neutral growth, 77 = T = 0 , as 
illustrated by the solid black line in Fig. [^ 

For 77 0 and T yf 0, species reproduce with different 

rates. Hence, the fronts of two neighboring sectors (occu¬ 
pied by different species) advance with different speeds. 
This influences the tilt of the front orientation, in ad¬ 
dition to stochastic roughening in neutral systems, and 
thereby the movement of the separating boundary. Thus, 
we do not expect Eq. § to be valid. 

Indeed, Fig. reveals different regimes for the mean 
number of boundaries Ni,. For T > Tc, Nj, decays expo¬ 
nentially in time in line with the exponential decay of the 
survival probability Pc in Fig.[^and similar to the case of 
selective advantage in Ref. [28] T For P < Tc boundaries 
annihilate less frequently. Narrow defector sectors persist 
in the front dominated by cooperators since almost all in¬ 
dividuals in the defector sectors have cooperating neigh¬ 
bors. This results in the upward curvature of the curves 
in Fig. However, the defector sectors cannot expand 
laterally and ultimately loose contact to the front due to 
random fluctuations, and Ni, declines exponentially. At 
Tc the number of boundaries decreases with a power law 
TVf, ~ where a new exponent y « 2.50±0.05 appears. 
This power law implies that the number of boundaries de¬ 
clines from the initial value Ni,{t = 0) ^ L to the order 
of 1 in the coarsening time 

rl/x ~ r0.40±0.01 
^coarse L/ L/ . 

Comparing with Eq. ([^ reveals l/y < This 

suggests that for large systems fixing the front to one 
species takes much longer than coarsening to a few sec- 















FIG. 9. Time evolution of two boundaries with initial dis¬ 
tance I //2 in a system close to criticality (T = 1 . 6 ) and at 
7?. = 0.1. Standard deviation a of the transverse distance plot¬ 
ted versus time t for several system sizes L. The boundaries 
move superdiffusively. Initially, a ~ with rj = 0.71 ± 0.05, 
and consistent with Eden scaling bold black line). At 

later times a crossover to cr ~ t’’ with rj' — 0.9 ± 0.1 occurs 
(dashed black line). 



FIG. 10. The width w of the front plotted versus time for 
several system sizes L close to the critical point at T = Tc ~ 
1.6. Inset: After rescaling time f by Tx ~ and width w 
by Wx ~ the curves collapse onto a single master curve. 

The solid black line indicates Eden roughening with w ~ 
whereas the dashed black line shows selective roughening with 
w ~ for t ^ Tx- Note that the saturation of w at large 
times is due to the finite system size. 


tors. Hence, the few remaining boundaries move differ¬ 
ently compared to early times since they have to annihi¬ 
late to fix the front to a single species. 

To check if this is the case, we employ the initial con¬ 
dition, where the front is composed of only two sectors of 
size L/2 each, separated by two boundaries. We quantify 
the boundaries’ random motion by monitoring the tem¬ 
poral evolution of the standard deviation for the trans¬ 
verse distance £±, 


a{t) := ■ (9) 

We subtract the mean distance {£j_{t)) to take care of any 
transient drift, when the front relaxes from its initially 
flat to the rough shape, and an expected small drift if T 
is not exactly Tc- 

From Fig. |^we see that, for early times, tr grows like 
a power law, a ~ with rj = 0.71 ± 0.05. This is 
consistent with meandering boundaries induced by Eden 
roughening, cr ~ . We expect such a behavior 

since the roughness of the front has not fully developed 
yet. For large systems and later times we find a crossover 
to cr ~ with rj' = 0.9 ± 0.1. This confirms our ear¬ 
lier statement that at late times the few boundaries re¬ 
maining after coarsening move differently. Indeed, they 
show an even stronger superdiffusive motion than Eden 
scaling, which is associated with the long-lived and pro¬ 
nounced surface undulations in systems with frequency- 
dependent selection. 


D. Surface roughening of the expanding front 

We now discuss the surface roughness or undulations of 
the expanding front and compare our results to the clas¬ 
sical Eden model. We measure the roughness of a front, 
which has not yet fixed to one species, by calculating its 
width w for system size L and at time t: 


/ , L _ ^1/2 

{L,t) := - hit)]"^ ] . 


( 10 ) 


Here, h(i, t) is the longitudinal position of the front at its 
transverse site i and h{t) is the mean position 


Ki) ■= ■ 


( 11 ) 


i=l 


Figure 10 plots the width w{L,t) close to the critical 
point at T = Tc ~ 1.6. Initially, the roughness of the 
front grows like in the original Eden model [62] : re ~ t'*', 
where the growth exponent 7 = 1/3 belongs to the KPZ- 
universality class [63] . At intermediate times a new 
regime sets in, where w V increases with enhanced 
growth exponent 7 ' « 1.3±0.1. This marks the transition 
from Eden roughening to “selective roughening”. Here, 
the typical shape of the front is determined by advancing 
cooperator sectors and trailing defector sectors and ulti¬ 
mately drives the accelerated increase in the front’s width 
w. The crossover to this regime happens at time Tx , 
which increases with system size L as Eig.|10|shows. This 
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makes sense since we expect selective roughening to dom¬ 
inate over Eden roughening when the lateral extension 
of sectors is comparable to L, i.e., Tx Tcoarse ~ 

Due to Eden roughening the width at the crossover is 
Wx ~ Tx ~ Indeed, rescaling width and 

time with Wx and Tx, respectively, collapses all data in 
Fig. onto a single master curve, as the inset demon¬ 
strates. To conclude, surface roughening close to critical¬ 
ity occurs in two regimes. Until crossover time Tx , one 
observes Eden roughening, whereas for times larger Tx 
selective roughening occurs until the front fixes to one 
species. The dynamics of the width of the front is sum¬ 
marized by 


w(E, t) 


t^Tx'- Ei/x , 
tijJa-i'Mx t-^Tx ■ 


( 12 ) 


V. SUMMARY AND CONCLUSION 

In this work we studied a generalized Eden model, 
where two species compete with each other at the rough 
expanding front. Individuals of the two species influence 
each other by frequency-dependent selection, which acts 
between nearest neighbors. We analyzed the evolution¬ 
ary dynamics at the expanding front, where single-species 
sectors form and coarsen. Ultimately, the front fixes to 
one species, which we identify with an absorbing state of 
our model. 

In its general form the model can implement several 
scenarios including selective advantage, and also well- 
known game theoretical settings like the snow drift game 
or the coordination game. Each of them creates distinct 
patterns, which should be analyzed in detail in future 
work. For the prominent example of a public goods 
game, we find that cooperators prevail in a wide param¬ 
eter regime, as expected for a spatial version of a social 
dilemma daiisss]. For other parameter values defec¬ 
tors take over the front, as usual. 

We identify the transition between long-term coopera¬ 
tion and long-term defection as a nonequilibrium critical 
phase transition between two absorbing states. The set 
of critical exponents (see Table |l|, which we determined 
by analyzing critical and finite size scaling, shows that 
the phase transition belongs to a new universality class. 
We attribute this result to the fact that the front in our 
model is rough and not flat as in usual absorbing states. 
Close to the critical transition the front’s roughness ex¬ 
hibits a crossover in time from slow Eden roughening 
to fast selective roughening. Strong roughening has also 
been observed at phase transitions in a related model by 
Lavrentovich and Nelson m- 

Interestingly, the critical exponents determined in our 
present work violate the so-called generalized hyperscal¬ 
ing relation m 


(13) 


where the number of transverse dimensions d is 1. This 
relation holds for most universality classes of phase tran¬ 
sitions to absorbing states m- It is not obvious why our 
model does not obey the scaling relation. 

The roughness of the front correlates with superdiffu- 
sive motion of the boundaries separating sectors. Two 
factors contribute to the movement of the boundaries on 
long length scales. On the one hand, the direct competi¬ 
tion between the species on either side of the boundaries 
pushes them towards the sector composed of the more 
slowly reproducing species. On the other hand, bound¬ 
aries follow the local tilt of the front. In the public goods 
game cooperator sectors are advanced, while defector sec¬ 
tors lag behind. Near the critical transition, defectors 
outcompete their direct cooperating neighbors but the 
front is tilted towards sectors filled by defectors, so the 
two factors move the boundaries in opposite directions. 
At the phase transition both effects cancel and the front 
fixes with equal probability to either species. The strong 
roughening correlates with superdiffusive motion of the 
boundaries with nearly ballistic scaling. 

Accordingly, whether a species takes over the expand¬ 
ing front is determined by two contributions: its repro¬ 
duction rate relative to its competitor and its position 
relative to the average front position. The influence of 
different reproduction rates of neighboring species can 
directly be compared and is summarized in the phrase 
“survival of the fittest”. The position at the front deter¬ 
mines the available space for progeny, which then have 
the opportunity to expand sidewards. This is illustrated 
by the phrase “survival of the fastest” [38]. 

In our model the number of sectors only decreases. It 
does not include experiments with mutually beneficial 
interactions between different species, which do not gen¬ 
erate sectors EZl IHOIIHI]. In future extensions of our 
model this may be remedied by including motility of in¬ 
dividuals (82] [83] , by allowing reproduction to more dis¬ 
tant lattice sites m, or by increasing the maximal num¬ 
ber of individuals per lattice site from one. Moreover, 
it is worthwhile to consider interactions ranging beyond 
nearest neighbors, since biomolecules, released by indi¬ 
vidual microorganisms, may diffuse in the extracellular 
medium In the public goods game scenario this 

would stabilize narrow sectors of defectors so that they 
do not lose contact to the front. 

In general, range expansion of multiple species will de¬ 
velop enhanced roughness at the growing front. As we 
demonstrated here, the corresponding models have new 
and interesting statistical properties. From a biological 
point of view, roughness is important. It affects the ter¬ 
ritories that different species occupy and thereby their 
evolutionary success through the strong random motion 
of sector boundaries. This may also be relevant for range 
expansion in a real environment and not just in a test 
tube. To better understand the properties and conse¬ 
quences of rough expanding fronts, further theoretical 
work is needed. At the same time further experiments 
should look for the fingerprint of roughness in microbial 
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